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Abstract 

Starting from relativistic quantum field theories, Kovtun et al. (2005) have 
quite recently proposed a lower bound rj/s > ft/(47rfcB), where i] is the shear 
viscosity and s the volume density of entropy for dense liquids. If their proposal 
can eventually be proved, then this would provide key theoretical underpinning 
to earlier semiempirical proposals on the relation between a transport coefficient 
rj and a thermodynamic quantity s. Here, we examine largely experimental data 
on some dense liquids, the insulators nitrogen, water, and ammonia, plus the 
alkali metals, where the shear viscosity rj(T) for the four heaviest alkalis is known 
to scale onto an 'almost universal' curve, following the work of Tankeshwar and 
March a decade ago. So far, all known results for both insulating and metallic 
dense liquids correctly exceed the lower bound prediction of Kovtun et al. 
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1. Background and outline 

Going back, at very least, to the early semiempirical proposal of Rosenfeld 
[1], quite an active field has grown up relating shear viscosity rj to the volume 
density of entropy s in dense liquids @, S S] • Some motivation for returning 
to this general area, however, has come from unexpected quarters, namely rela- 
tivistic quantum field string theory methods. Thus, Kovtun et al. [f|, starting 
from such an ah initio viewpoint, have proposed that the ratio rj/s has a lower 
bound given by fi/(47rfc B ) ~ 6.08 • 10~ 13 K-s. Here, h and fc B denote Planck's 
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and Boltzmann's constants, respectively. While the lower bound proposal is 
not yet proven, their arguments are highly plausible (see Section [5] for a brief 
summary) . 

The outline of the present paper is then as follows. After a brief review 
of the hard sphere model in Section [21 with specific emphasis to its results 
for the shear viscosity and entropy density in dense fluids, in Section [3J we 
examine some scaling properties of the three dense insulating liquids He, N2, 
and H 2 considered in ^j, as well as NH 3 . This is then followed in Section 0] 
by work on dense liquid metals [f|. Some six liquid transition metals treated 
computationally by means of the embedded atom model are then considered, 
with the alkalis playing the dominant role because of already known scaling 
properties of the transport coefficient rj(T) for the four heaviest alkali metals, 
Na through Cs Q , these results constituting Section 14.21 A discussion follows 
in Section [5l together with proposals for further work that should prove fruitful. 



2. Hard sphere model of shear viscosity rj and entropy density s in a 
dense fluid used in phenomenological context 

We begin with what we believe is the simplest model of a dense fluid, viz. 
that composed of hard spheres (HS). In early work, Collins and Raffel [|| and 
Longuet-Higgins and Pople |9| stressed that in a liquid of rigid molecules the 
singular nature of the interaction allows a finite flux of momentum and energy 
even when the liquid radial distribution function is momentarily isotropic, as it 
is in equilibrium. 

The results obtained by Longuet-Higgins and Pople can be displayed in 
terms of the deviation of the equation of state from the ideal gas form as 



?7hs = tP g 
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P 
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where a denotes the hard sphere diameter, while M is the particle mass. The 
factor in square brackets in Eq. ([1]) embodies the probability of finding two 
spheres at contact, which in turn is given by the contact value of the pair 
distribution function g(r). 



The Percus-Yevick theory [lfj, [11|, [12| of the hard sphere fluid leads to a 
number of analytic, though of course approximate, results, which will be used 
below in exemplifying the hard sphere dense fluid predictions pertaining to the 
Kovtun et al. [5| lower bound on the ratio shear viscosity 77 to volume density 
of entropy s. 

The so-called excess entropy and other thermodynamic properties of the 
density hard sphere fluid can be usually derived with adequate accuracy from 
the Carnahan-Starling [l3| equation of state. In terms of the packing fraction 
Pf, this reads 

PN k B r (i- Pf r ' 
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where pf — irpNcr 3 /6, with pN = N/V being the number density. Eq ^ can 
be generalized to embrace hard fluids of non-spherical molecules as [14| 

P 1 + (3a - 2)p f + (a 2 + a - l)pj - a(5a - 4)p 3 f 



p N k B T (l~p f ) 3 

with a = 2.2346 for tetrahedral molecules, and a = 1 for spherical molecules, in 
which case it reduces to Eq. ([2]) . Hence, by specifying the pressure P as in the 
examples considered by Kovtun et al. [5(, number density pn and temperature 
T are related via the packing fraction pf by Eq. ([2|). Insertion of Eq. |2]) in the 
right-hand side of Eq. |T]) then yields t^hs — Vasi 171 ! T; a). 

To construct the Kovtun ratio r\j s in this model, we next write the excess 
entropy Se in the form 

Srs _ Se — Sp _ 2 1 3 — — (4) 

k B k B l — Pf (1-J 3 /) 2 fe' 



where So is the entropy for an ideal gas [15j . 
So 5 
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const + log T 3/2 - log p N , (5) 



the latter functional form, including an additive constant, being derived through 
thermodynamic considerations only e.g. by March and Tosi [1 61 ] . 

The excess entropy Se in Eq. ((4]) is clearly positive in the range < pf < 1/2. 
Dense fluids, for orientation, freeze when pj ~ 0.46. The volume density of 
entropy shs — pSHs/k B for the dense hard sphere fluid is hence known from 
Eq. (j4}. The ratio rj/s can then be formed from Eqs. JT]) and (j4]), and plotted as 
a function of T, after eliminating pf through the hard sphere equation of state 
Eq. @ choosing the pressure P — 100 MPa [f| (Fig. [I}. Appendix [K\ gives a 
formal proof relating to r\j s under isobaric conditions, that 77/s is a function of 
T, as in Fig. [T] (see also Fig. HJ). 



3. Shape similarity after scaling of available experimental data on 
dense insulating liquids helium, nitrogen, water, and ammonia 

Though our main focus will be on dense metallic liquids, it seemed natural 
to start by examining scaling properties for helium, nitrogen, and water, the 
insulating liquids already considered by Kovtun et al. [5| along with their pi- 
oneering proposal. In their Fig. 2, all three liquids are treated in their plot of 
the ratio rj/s vs absolute temperature T. 

Because of the lower bound proposal, we have replotted their data in Fig.[2k. 
This scales T with T min , and also scales (t]/s)t in units of the value at the 
minimum, (?7/s)x mln • It is of interest that for these very different liquids, the data 
collapses somewhat to exhibit considerable shape similarity. For completeness, 
as well as to gain insight as to the source of the minima in Fig. [2^,, we show 
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Figure 1: Schematic scaled data for ratio rj/s for P = 100 MPa, within the hard sphere model. 
Data sources referred to in 

s(T) and rj(T) separately in Fig. [2Jd for nitrogen, water, and ammonia. Over 
the range of temperatures displayed there, the viscosity is featureless, and it is 
the entropy density s(T) that has the characteristic features. It is noteworthy 
that the maximum in s(T) for water is near to the critical temperature T c . This 
is also true for nitrogen. 

Having examined these four dense insulating liquids, we turn to our central 
examples, an embedded atom model computer simulation of six liquid transition 
metals, and, following that, the dense metallic alkali liquids Li through Cs. 



4. Some relevant results concerning transport in liquid metals 

So far, the specific dense liquids for which we have presented experimental 
results are insulators: atomic He and molecular fluids H2O and N2. Though we 
do not have such complete experimental data, we feel it of interest to summarize 
here some obviously relevant results on some twelve or so liquids having itinerant 
electrons, viz. liquid transition metals (Cu, . . . ) and the five liquid alkali metals 
Li, Na, K, Rb, and Cs. 

4-1- Computer simulation results on some liquid transition metals 

Hoyt et al. have carried out important computer simulation results of 
liquid transition metals, using as interatomic potential for their studies those 
generated by the embedded atom model (EAM). The transport coefficient they 
focussed on was the self diffusion D. They were able to refine, for Cu, Ag, Au, 
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Figure 2: (a) Scaled data for ratio r)/s for insulating dense liquids He, N2, H2O, and NH3. 
Data sources referred to in The resulting values (in SI units) for (??/s) m i n are 25.07 for 
He, 24.54 for N2, 45.15 for H2O, and 94.00 for NH3. (b) Separate variations with temperature 
T of volume density of entropy s and of viscosity r\ (right-hand scale) for nitrogen, water, and 
ammonia. Data sources in . (c) Correlation between critical temperatures and temperatures 
T m i„, at which rj/s is minimum in Fig.[2ji (closed symbols), and temperatures Tmax, at which 
the entropy density is maximum in Fig. [2p (open symbols), for nitrogen, water, and ammonia. 
The straight line is a guide to the eye. 



Ni, Pd, and Pt, the earlier important study of Dzugotov 2], who approximated 
the excess entropy via the pair distribution function g{r). Hoyt et al. |3j found 
a form like that of the semiempirical proposal of Rosenfeld [17J , who wrote (see 
also Ref. 0) 

D ~ 0.6/9 w 1/3 wexp(-0.8s), (6) 

where s = —SE/Nk^ 1 with Se the excess entropy. Hoyt et al. Q recovered 
a form like Eq. © for the six transition metals listed above, but with the 
factor cxp(— s) to good accuracy. The scatter of Dzugotov's results B was 
substantially reduced by the computer simulation results of Hoyt et al. [3| . 

As suitable results for excess entropy are not available to us, we next note two 
points: (i) The Stokes-Einstein relation (see Ref. can be used to estimate 
the viscosity of the above transition metals from the Hoyt et al. form for D, 
and (ii) it is of obvious interest to examine the relation of the above to the 
predictions of the hard sphere model discussed above. This is because of the 
early study of Bernasconi and March [l8| (see also Ref. These authors 

calculated from experiment near freezing the direct correlation function c(r) at 
r = 0, which for hard spheres the Percus-Yevick theory invoked above predicted 
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Figure 3: Experimental results for the number density pjy as a function of temperature T 
for H2O and NH3 under isobaric conditions (P = 100 MPa). Data sources referred to in 
Behavior of modified hard sphere result in Eq. J3} for a = 1 and 1.7, with a and hard sphere 
diameter a = 1 — 5 A treated as parameters, is also shown. 



Metal 



-c(r = 0) -c(k = 0) 



Cu 
Ag 
Au 
Ni 



GO 
51 
35 
41 



47 
53 
38 
■50 



c(r = 0) 
c(k = 0) 



1.3 
1.0 
0.9 

0.8 



Table 1: Direct correlation function c(r ; 
transition metals. After Ref. Hal. 



: 0) and its Fourier transform c(k = 0) for four liquid 



to be related to its Fourier transform, say c(k = 0), by [20j 

c(k = 0) = 1 + c(r = 0). (7) 

Table [1] shows the values obtained by Bernasconi and March p[ for four of 
the six liquid transition metals considered in the EAM studies of Hoyt et al. [|| . 

A necessary condition for a liquid metal to be 'hard-sphere'-like is for the 
ratio c(r = Q)/c(k — 0) to be near unity, as stressed by Bernasconi and March 
[]J]. This is well satisfied for Ag, Au, and Ni (Tab. [J). The polyvalent liquid 
metals Ga, Pb, and Sb are obvious exceptions from the study of Bernasconi and 
March [Tsj ]. the ratios being respectively 0.2, 0.4, and 0.3. 

Because of the results in Table [TJ we have been motivated to return to 
the hard sphere model. However, we have transcended the analytic model of 
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Longuet-Higgins and Pople [§] used above, by invoking computer simulation 
results. Thus, Speedy 2l| has fitted molecular dynamics results for D in the 



hard sphere model by the approximate formula 

D = D {l-^)[l + x 4 (0A-0Mx% (8) 

where x = pno~ 3 is proportional to the packing fraction pf introduced above, 
while 

3a fk B T^ 1/2 



D = -p- — (9) 
is the infinite dilution value of D, which is known exactly from kinetic theory. 



For comparison with Eq. © due to Speedy [211] . an alternative expression 



for the density dependence of the self-diffusion for hard spheres has been given 



by Erpenbeck and Wood |22| as 

D = D E (1 + 0.0382x + 3.18x 2 - 3.869a; 3 ), (10) 

where 

D E = i21^, (11) 

9hs(p + ) denoting the hard sphere pair distribution function at contact, referred 
to earlier in relation to the Longuet-Higgins and Pople study. This is related 
to the packing fraction pf by 

*-<- + > = ^- < 12 > 

Eqs. ([5]) and pH)) arc then compared in Fig. 2] as a function of packing fraction 
Pf. Since our main focus here is the shear viscosity 77, we have used the, of 
course approximate, proposal of Longuet-Higgins and Pople [9] for the specific 
form of the Stokes-Einstein relation, viz. 

Drj = -^<j 2 PN k B T, (13) 

to plot 77 vs pf in Fig. [5] 

4-2. Scaled experimental data on liquid alkali metals 

To conclude this section, it seemed of interest to refer to scaling properties 
of the five liquid alkali metals Li, Na, K, Rb, and Cs 0- Why the liquid alkali 
metals are especially important in the present context is clear from the scaling 
exhibited in [7j of the shear viscosity r](T) from available experiments. In fact 
the so-called fluidity v^ 1 is the quantity plotted, but appropriately scaled by its 
value at the melting temperature T me i t , to distinguish from the minimum T m i n 
in the plot of the ratio n/s vs T (Fig. 2]) 

f.0i 



We have redrawn Fig. 3 of Ref. 



in Fig. [6] and it is quite clear that, for 



the four heavier alkalis, Na, K, Rb, and Cs, the scaling used collapses the data 
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Figure 4: Normalized self-diffusion coefficient D/Do as a function of packing fraction py within 
the hard sphere model, according to the Speedy [Eq. ©] and Erpenbeck and Wood [Eq. JlOt l 
results. 



onto an 'almost universal' curve. While there is some shape similarity for the 
remaining (lightest) alkali metal Li, this has an appreciably different magnitude, 
under scaling, from the other four. In principle, one can estimate the volume 
density of entropy s as a function of T from experimental structure factor data. 

In fact, in a further paper [171 ]. Rosenfeld writes an empirical formula relating 
shear viscosity rj and scaled entropy sr, where 

s R = -(S-S )/NkB, (14) 
So referring to the ideal gas as in Eq. ©. Writing 

Vred = VPN V3 /(Mk B T)^ 2 , (15) 

where pat is the particle number density, while M is the ionic mass, Rosenfeld 
proposes the semiempirical formula 

r/ red = 0.2exp(0.8s R ). (16) 

Hence, it follows [13] that 

^=0.2^ 3 (Mfc B T)V 2 eXp( ° g 8SR) . (17) 

We have employed the experimental viscosity data in Ref. [2J] for liquid Li (see 
also Fig. [6} , to extract the Rosenfeld scaled excess entropy sr from Eq. (jTT]) , 
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Figure 5: Approximation to hard sphere viscosity r/HS using computer results for D/Do in 
Fig.g] plus Eq. Ifl3i 



when combined with the approximate equation of state of Chapman and March 
[23|. This reads 

PnKbT Pn^bT 

where T is the coupling parameter of the so-called one component plasma (see 



25() defined by 

r =^T' < 19 > 

where p^ 1 = 47ra 3 /3. Fig. [7)3 shows schematic plots of vs T, which it would, 
of course, be of considerable interest to have quantitative experimental results 
in the future at the thermodynamic conditions in Fig. [5] and in the experimental 
study of ij(T) in Ban et al. |24[. Fig. [7): finally shows the schematic form of the 
volume density of entropy s(T) obtained from Eqs. ([6]) and ([5]). 



5. Discussion and future directions 

Motivated by the proposal of Kovtun et al. 0] of a lower bound for the ratio 
of shear viscosity r\ to volume density of entropy s, we have studied this ratio for 
some fifteen dense liquids. These include liquid metals, five liquid alkalis, plus 
some six liquid transition metals, for the former of which the shear viscosity 
rj{T) has been determined experimentally over a range of temperature T. Using 
the scaling of these data in [7|] , Figure [6] has been constructed by making use of 
the semiempirical Eq. (fT6|) given by Rosenfeld [3, [I3]- Into this equation, the 
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Figure 6: Scaled inverse of shear viscosity r}(T) as a function of temperature T for liquid alkali 
metals (redrawn from Ref. 0). Note that, apart from the lightest metal Li, the other data 
fall on an 'almost universal' curve. Note that T m is the melting temperature. 



experimental data for rj(T) has first been used to predict the volume density 
of entropy s(T) plotted in Fig. [Jj It will be of significance for the future to 
refine these data using X-ray or neutron diffraction measurements to estimate 
semiempirically the pair correlation function g(r). This is then sufficient to 
determine the entropy s, following, for instance, the procedure of Dzugotov 

The most quantitative results in the present study are presented in Fig. O 
In part (a), we have scaled data for the four insulating liquids shown by plotting 
the ratio rj/s vs T/T m { n , where T m j n is the temperature at the minimum of 77/s 
in each case. It is seen that there is considerable shape similarity for the four 
liquids, in spite of the range of isobaric conditions involved. Fig. [lb is presented, 
again from experiment, to give insight into the origin of the minima in rj/s. Over 
the range plotted in Fig. v(T) is monotonically decreasing with increasing 
temperature, and the minimum in rj/s therefore arises from the turning point 
in the volume density of entropy s(T), at temperature T max . Fig. [2}: shows that 
both T m i n and T max are of the order of the critical temperature T c for N2, H2O, 
and NH 3 . 

The remainder of the article begins a study of the ten liquid metals already 
mentioned. Because of the incomplete availability of (i) experiment and (ii) 
computer simulation data, we have had recourse to models. These are (i) the 
well known HS model, set out in Section [21 where both approximate analytic 
theory and computer simulation data is now available (see especially Fig. [4] for 
the diffusion coefficient D). This has been combined with the (of course ap- 
proximate) analytic form of the Stokes-Einstein Eq. (|13p to construct Fig. [5] It 
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would, of course, be valuable to have results for the shear viscosity rj paralleling 
those displayed in Fig. 0] for D, to test the HS predictions of Fig. [5l and if 
necessary to refine them. 

We may mention also at this point the gen eralization Eq. © of the HS 
equation of state ([2]). It has been suggested 14 1 that Eq. ([3]) may have some 
relevance to water with a ~ 2, which is close to the numerical value 1.7 employed 
to construct Fig. [31 relating to experiments on water and ammonia. While on 
the HS model, the 'softening' of the core to inverse power repulsive potentials 
is referred in Appendix [Bl (see also especially Fig. [5]). 

For the ten or so liquid metals considered, and beginning with the six tran- 
sition elements studied by Hoyt et al. [3| using the embedded atom model, 
Table [1] shows that a necessary, though not sufficient, condition for HS behavior 
(that the final column is equal to unity) is well obeyed for two, Ag and Au, and 
is approximately valid for Cu and Ni. So in the future, it will be interesting to 
bring computer simulation and experiment into contact with the HS model, and 
especially with the predictions for the shear viscosity r\ in Fig. [5] Finally, for the 
five liquid alkalis, Fig. [6] summarizes scaling properties of the shear viscosity. 
However, experimental results on the density-temperature relation under the 
thermodynamic conditions appropriate to Fig. [5] and to the isotope experiments 
of Ban et al. [24] would be most helpful in the present context, as, next to of 
course He (see Fig. [2J), Li is the most likely candidate for quantal effects, and 
thus for a possible lower bound to compare with the prediction of Kovtun et al. 
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In a related vein, we note that the work of Tosi [26[ on electron plasma 
theory has led to the conclusion that the itinerant electrons in the alkali metals 
contribute less than 10% to the shear viscosity. As a proposal for future work, 
though somewhat esoteric, it would be of interest to study in its own right the 
homogeneous electron liquid [27j |. 

As for liquid Li, we believe it would be of considerable interest to use a 
density-dependent pair potential, such as that derived by Parrot and March 
[28j from first principles electron theory, to calculate the excess entropy ap- 
proximately from the resulting pair distribution function g(r), as set out in [2|. 
Then, adding the ideal entropy in Eq. ([5]), one could plot the volume density 
of entropy, say, along an isobar. This could then be compared for any shape 
similarity with the three insulating liquids treated in Fig. [2}d. 

As this study was nearing completion, we became aware of the highly relevant 
work of Schafer [ijj ]. This author studied the ratio rj/s for trapped fermions in 
the unitarity limit. The result of this investigation was the prediction that 
rj/s is roughly 1/2 in units of h/hs. This example prompts us to conclude this 
article by summarizing briefly the uncertainty principle argument pointed out by 
Kovtun et al. [f|. This consists of three parts: (i) r\ for a plasma is proportional 
to er, where e is the energy density while r is the typical mean free time of 
a quasiparticle, and (ii) the entropy density is proportional to the density of 
quasiparticles, say n, i.e. s ~ fceii. Thus the desired ratio rj/s ~ ks~ er/n. 
The third part of the argument [B[ invokes, as mentioned above, the uncertainty 
principle in the form (e/h)r > H, without which the concept of quasiparticle 
lacks meaning. This gives rj/s > H/kn, without, however, the factor 1/Att in ||. 
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A. Scaled diffusion coefficient in terms of excess entropy Se in the 
hard sphere model: Relation to the work of Kovtun et al. 



In Fig. [5] we compare the computer simulation results of Speedy 21[ with 



those of Erpenbeck and Wood [22j. Evidently, with refinements to available 
computer data, the two curves in Fig. [4] will converge to a single curve which 
we represent by writing the scaled diffusion coefficient D/Dq in the functional 
form 

£ = /(*>/)• (20) 

But for this HS model, we know the excess entropy Se in the Percus-Yevick 
approximation as is given in Eq. From Fig.[3]representing an approximation 
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to the formally exact Eq. (j2"0)) , we can conclude that a unique inversion exists, 
to yield 

Pf = g(D/D ). (21) 

Inserting Eq. (|2ip into Eq. we therefore conclude that there exists the 
functional relation 

= h(D/D ) (22) 

in the HS model. 

If we combine this exact formal result, Eq. (|22[) with the (now approximate) 
Stokes-Einstein relation proposed by Longuet-Higgins and Pople [9(] for the HS 
model, Eq. (|13[) . and use pf = irpNcr 3 /6 to find 

Dr,= ^k B T, (23) 

we find from Eqs. (gDJ) and (HU that 

But with the plausible assumption that Eq. (|22p has the unique inversion 

° >m (25) 



Do V * 

wc find by combining this equation with Eq. (|24p that 



wA^AV, (2g) 



We propose that Eq. (f2T))) demonstrates the functional form of in terms of the 
dimcnsionless excess entropy Ss/ks in the HS model. 

But when we turn to consider the ratio r\j ' s as in the lower bound proposed by 
Kovtun et al. [fa], we note that the total entropy S is given by Eq. (fl4|) . While 
SW&b in Eq. ([22]) is solely a function of the packing factor pf, Sq involves 
also, following Landau and Lifshitz [15|, the thermal de Broglie wavelength A = 
(2irh 2 / Mk^r) 1 / 2 , providing a natural length unit for the volume per particle 
pjf 1 , and hence the volume entropy s = p^S depends on a, pf, and A, as does 
the free energy F since the ideal gas result Fq reads [IH 

F = -Nk B T[l - log^A 3 )] (27) 

In the isobaric plots of Fig.0 pn = Pn(T), and one has that -q/s = n(T). 

B. Results for inverse power repulsive potentials 

In the main text, together with experiments and their interpretation espe- 
cially on the molecular insulating liquids nitrogen, water, and ammonia, we have 
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Figure 8: Entropy for water and nitrogen at the given pressures (data source as in Kovtun et 
al. §), as a function of p/T 3 / n , for n = 12. Data for water have been shifted vertically for 
clarity. 



given prominence to the hard sphere model. Here, we generalize this model to 
the case of inverse power repulsive potentials having the form 



(r)=e^) n . (28) 



The merit of this model is that, for any such inverse power potential, thermody- 
namic properties are readily accessible since only a single isotherm, isochoric, or 
isobar needs to be known. Then all other properties can be found, as discussed 
by Hoover and Ross [13] ■ 

Quite specifically, for pair- wise additive potentials having the form (|28|) . 
(now dimensionless) thermodynamic properties depend only on the one density- 
demperature variable, y say: 

v = p (lkT> (29) 

where p is proportional to the packing fraction pf and is defined as p = Na 3 /V. 

We have studied the data on nitrogen and water presented in Sec. |3]of the 
main text, and for n = 12 we find remarkable shape similarity (Fig. [8|). 

Because of transport coefficients being another focal point here, we finish 
this Appendix by referring to the interesting study of the diffusion coefficient 



D(n) for the potential (|28|) by Heyes and Powles 3l|. For a packing fraction of 



0.044, these authors studied how D{n) approached the hard sphere limit D(oo) 
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and found that their data could be well fitted by the form, with D now in 
dimensionless form, 

D(n) -D(oo) = const ■ n -1 " 71 , (30) 
indicating a rather slow approach to the hard sphere result. 
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